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ABSTRACT 



We examine a simulation of flux emergence and cancellation, which shows a complex 
sequence of processes that accumulate free magnetic energy in the solar corona essential 
for the eruptive events such as coronal mass ejections (CMEs), filament eruptions and 
flares. The flow velocity at the surface and in the corona shows a consistent shearing 
pattern along the polarity inversion line (PIL), which together with the rotation of the 
magnetic polarities, builds up the magnetic shear. Tether-cutting reconnection above 
the PIL then produces longer sheared magnetic field lines that extend higher into the 
corona, where a sigmoidal structure forms. Most significantly, reconnection and upward 
energy-flux transfer are found to occur even as magnetic flux is submerging and appears 
to cancel at the photosphere. A comparison of the simulated coronal field with the 
corresponding coronal potential field graphically shows the development of nonpotential 
fields during the emergence of the magnetic flux and formation of sunspots. 

Subject headings: MHD — Sun: interior — Sun: atmosphere 



1. Introduction 

Regions of intense solar magnetic fields such as sunspots and active regions are known to exhibit 
ener getic outbursts that are manifest in many forms, such as CMEs , filament erupti ons, and flares 
(e.g. iForbesj hood : IlowI boOll : iForbes et all bood : iGopalswam^l bood : ISchriiverl bood ) . In the most 
dramatic form, magnetic fields release energy up to 10 32 ergs, and in the process ex pels 10 15 gm of 
plasma and 10 20 Mx of ma gnetic flux into interplanetary space (IBieber &; Rustlll995T) at speeds that 
can reach 3000 km/s (e.g. iGosling et al1li99d : 1st. Cvr et aDbood : Izhang et allbood ). Although 
less energetic, large-scale regions of relatively weak magnetic fields can also erupt giving rise to 
streamer blowouts and the expulsion of quiescent filaments into the solar wind (e.g. 
199.4 *~ 



Hundhausen 



Lvnch et al.ll2010h . 



In order to produce such eruptions, the magnetic field must possess free magnetic energy, which 
requires that the field is in a nonpotential configuration with electric currents passing through the 
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photosphere and into the corona. In general, the field lines of the non-potential force-free fields 
are stretched, expanded and intertwined compared to the equivalent potential field mapping to the 
photosphere. Such force-free fields in the solar corona have historically been taken to be of two 
forms: arcades, in which the foot points have been sheared to produce a strong field component 
parallel to the polarity inversion line (PIL), or flux ropes, in which the magnetic field is twisted 
about a central axis. This dichotomy of field configurations naturally leads to two broad groups of 
theoretical and numerical models of CMEs. 

In the case of magnetic arcades, initial configurations are typically taken to be line-tied to the 
lower boundary and subsequ ently energized with prescribed shear flows, which cause the arcade to 
expa nd and then erupt (e.g. ISteinolfsonlll99ll ; iMikic fc Linker! 1 19941 ; lAmari et al.lll996l ; iGuo Sz Wu 
1998). Another variant, the so-called breakout model, involves a system of quadrupolar fields 
with a system of three arcades. The field reverses direction over the central arcade allowing mag- 
netic reconnection to rele ase the central arcade (e.g. Antiochos et al.lll999l ; Lynch et al. 20041 . 12OO8I ; 
van der Hoist et al.ll2009l ). In all of these cases, magnetic flux ropes form during the eruption pro- 
cess and are expelled from the corona. For those models, which assume that flux ropes exist prior 
to the eruption, the flux ropes are assumed to reside within an arcade where the outward-directed 
Lore ntz force of the rope (hoop force) is balanced by magnetic t ension of the surrounding arcade 



(e.g. iForbes h Isenberg 



1991 



Low 



1994 



Titov Demoulinl ll999). The central question with flux 



rope models is how the force imbalance occurs and results in eruption. Fl ux ropes have been ener- 



Tokman fc Bcllan 



2002 



gized by the applicatio n of rotational motions at the rope footpoints (e.g 
Rachmeler et al.ll2009i ) that twist up the field. Flux ropes have also been forced to emerge into the 
corona by bodily adyecting the m through the lo wer boundar y of th e computational domain (e.g. 
Torok & Klie ml I2OO5I : iFanlbood ). In the case of kmari et al.l ((2003), flux cancellation at the PIL 
forms a flux rope within a sheared magnetic arcade. Following these formation processes, the flux 
rope may lose equilibrium and erupt from several in stabilities including: unbalance of magnetic 



forces, which is foll owed by magnetic reconnection (Forbes &: Isenberg 



1991; Amari et al 



2003 



Rousscv ct al. 


2004 


), torus instability ( 


Torok & Kliem 2005 


), or kink instability ( 


Sturrock et al. 


2001 




Fan 


200 


A. 







Plentiful observations provide evidence for both arcades and flux ropes as the coronal field 
structures for CMEs. S heared magnetic fields are measured directly at the photosphere with vec- 
tor magnetogra phs (e.g. lHagyard et al.lll984l : IZirin k. Wanglll993l : lFalconerll200ll : lYang et al.ll2004l : 
Liu et al.ll2005r). where t he field is found to run nearly parallel to the PIL. More recent analysis 
by ISchrijver et al.l (|2005l ) found that shear flows associated with flux emergence drove enhanced 
flaring. Simil arly, active reg i on CME pro ductivity is also strongly correlated with magnetic shear, 
as shown by IFalconer et al.l ([2002], 120061 ). There is also evidence suggesting that flux ropes may 
bodily emerge through the photosphere to reside in the corona. The twist of the emerging flux 
tube is characterized by the presence of elong at ed magnetic po l aritie s, i.e., magnetic tongues in 



longitudinal magnetograms (ILuoni et al. 



2011 



Leka fc Steinerl (|200ll ) examined a time series of 



vector magnetograms of emerging magnetic flux and found that magnetic fields passing through 
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the photosphere are already twisted. Similarly, iLites et al.l (|1995l ) examined a time series of photo- 
spheric vector magnetogram s, whi c h the y found could be fitted by the emergence of a spherically 
shaped magnetic flux rope. ILites! (|2005l ) shows an example of the photospheric vector magnetic 
field and chromospheric structure seen in the Ha line, which suggest the presence of a coronal flux 
rope associated with a prominence. 

Magnetohydrodynamic (MHD) numerical simulations may provide new insights to our under- 
standing of the physical processes that build up the free magnetic energy necessary for solar erup- 
tions. Global scale simulations address how magnetic fields generated at the tachocline pass through 
the convection zone and have explained and reproduced many large-scale aspects of sunspots (e.g. 
Spruit fc van Ballegooijenl Il982l ; lAbbett et al.l l200ll ; iFanl 120081 ) . However because of the anelastic 
approximation, these models can only extend up to a height of 20 Mm below the photosphere. 
Above this depth, fully compressible MHD models have been employed to simulate the emergence 
of magnetic fields from the near surface convection zone into the lower corona. Earl y work demon 
strated the basic processes by which fields emerge to form bipolar active regi ons (e.g. IShibata et al 

These simulations were followed by cases (e.g. Manchester! 



198 



9: 



Matsumoto et al 



2001: 



2001 



Fan 



Magara &; Longcopel 2003 ). which showed that the Lorentz force, arising when fields expand 



into the highly stratified solar atmosphere, drives shear flows along the PIL. The shear flows then 
transport m agnetic energy into the c oronal portion of the flux rope which then expands, and drive 
its eruption ([Manchester et al.l 120041') . T hi s shearing process is suggested as the energy sour ce for 
CMEs and flares dMarichesterl l2007l. 120081) . iGalsgaard etaO (I2OO7I ) . I Archontis fc Torokl ||2008l ). and 
MacTaggart fc Hoodl (|2o"olF pr oduce fast CME-like eruptions by also invoking magnetic reconnec- 
tion with emerging flux ropes. Fan (j2008l ) simulated flux rope emergence that exhibited sunspot 
rotation driven by a tor sional form of the Lorentz force, which twisted up the coronal portion of 
the field as predicted by iParkerl (119771 ). However, in this case, no eruptive behavior was found. 



Recently, computing power has reached a level that has allowed the development of realistic 
solar models, including radiative and ther modynamic processes nec es sary to s imula te convection in 
conjunction with the upper atmosphere. Stein Nordlund ( 2006 ). Abbett (j2007l ). Rempel et al 



(|2009l ) , ICheung et al.l (|2010l ) , iKitiashvili et al.l (120101 ) , and lRempell (1201 if) emphasize the i mpor tance 
of the i n terac tion between magnetic fields and convective motions. IStein &: N ordlund (2006 ) and 
Abbettl (120071 ). have each addressed quiet-Sun magneto-convection, while IKitiashvili et al.l (120101 ) 
report the critical role of strong vortical d owndrafts around sm all m agnetic structures in the forma- 
tion of large-scale structures. The work of iRempel et al.l (|2009i ) an dlCheung et al.l (120101) t reats the 
formation and evolution of sunspots, including the Evershed flow. iFang et al.l (|2010l . 120121 ) address 
the emergence of magnetic flux ropes from a turbulent convection zone into the corona and found 
both shearing and rotational flows driven by the Lorentz force. These horizontal flows were found to 
dominate the ene rgy transport from the convection zone into the corona. Of particular significance, 
the simulation of IFang et al.l (120121 ) exhib ited a case of large-scale flu x cancellation, a phenomena 



strongly associated with CME initiation (jSubramanian &: Derdl200ll ). Here we examine the flux 



cancellation in conjunction with energy transport from the convection zone to the corona. 
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The paper is organized in the following way: Section [2] describes the numerical simulation; 
Section [3] studies the build up of magnetic shear in a flux cancellation event; followed by analysis 
on the free energy in the corona in Section 21 The energy transfer at the photosphere is discussed 
in Section [5j Section [6] summarizes our conclusions and discussion. 



Numerical Simulation 



To simulate the flux emergence in the convection zone, we first generate a relaxed solar atmo- 
sphere of dimension 30 x 30 x 42 Mm 3 , with the photosphere located at Z = Mm. Our model 
solves the MHD equations in Block-Adaptive Tree Solar-wind Roe Upwind Scheme (BAT-S-RUS) 
with additional energy source terms to describe thermodynamic processes in the sola r atmosphere , 



i.e., radiative co oling at the photosphere and in the corona, and coronal heating (jAbbett 



2007; 



Fang et al.ll2010l ). Meanwhile, implementation of the non-ideal tabular equation-of-state ([Rogers 
2000) provides a more accurate description of the partially ionized plasma in the convection zone. 
We use periodic horizontal, closed upper boundary conditions, and fix the density and temperature 
values at the lower boundary while setting the vertical momentum to be reflective. 

Taking advantage of the adaptive grid in BAT-S-RUS, our model produces a relaxed convection 
zone with a depth of 20 Mm, with a overlying photosphere of 0.45 Mm thickness and a corona of 20 
Mm, in a Cartesian do main. The v e rtical statification of the atmosphere in our simulation domain 
is shown in Figure 1 in iFang et al.l (|2012l ). We then carry out a simulation of the emergence of a 
buoyant, in i tially stationary, horizontal flux rope inserted at Z = — 10 Mm, shown in Figure 2 in 
Fang et al.l (|2012l ). Interaction of the rising flux rope with large-scale convective motion produces 
the bipolar structure of the flux rope, with the convective downflows fixing t he two ends of th e 
emerged section of the flux rope in the convection zone, illustrated in Figure 5 in lFang et al.l (|2012l ). 
Near-surface small-scale convection produces the magnetic polarities by coalescence and intensifies 
th e strength of the m agnetic flux. A small active region forms on the photosphere, seen in Figure 6 
of IFang et al.l (|2012l ). exhibiting strong interaction of the magnetic flux and the surface flows, i.e., 
the shearing and converging motions, separation and rotation of the magnetic polarities. In the 
pr esence of these flo ws, a flux cancellation event takes place at time t = 05:00:00, shown by Figure 6 
of lFang et al.l (|2012l ). during which 10% of the total photospheric unsigned flux is cancelled. During 
the flux emergence, horizontal motions dominate the energy transfer to the corona, while vertical 
flows transport energ y back into t he co nvection zone. Details on the dynamics of the emerging flux 
rope are discussed in IFang et al.l (|2012l ). 



3. A case of Flux Cancellation 



Here, we examine a region of magnetic flux cancellation at the photosphere in our simulation, 
which is also associated by high magnetic shear. Flux cancellation of this form is particularly 
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Fig. 1. — : The photospheric structure of B 2 in the area of flux cancellation at time t = 5:00:00 (a), 
5:10:00 (b), 5:22:00 (c). Blue and red arrows show the horizontal velocity in positive and negative 
magnetic polarities, respectively. Shear flows are most apparent at time t = 5:10:00. 
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Fig. 2. — : Left: temporal evolution of the total (blue), potential (red) and free (black) magnetic 
energy. Right: vertical distribution of free energy density at time t = 05:05:00. The dashed line 
indicates the altitude up to where 50% of the coronal free energy is stored. 



significant as it has been o bserved to be associated with flares ([Martin et al.l Il985l ) and CMEs 
(jSubramanian &: Derdl200ll ) in active regions. A clear example is observed in AR 10977, where re- 
connection at the PIL during flux cancell ation forms a highl y sheared arcade and sigmoidal structure 



that subsequently reforms after eruption (jGreen et al.ll201ll ) . iLi et al.l (120041 ) found converging flows 



consistent with flux cancellation in decaying and young active regions producing CMEs. Inspired 
by similar observations, converging motion and flux c ancellation have been imposed as bo undary 
conditions in MHD simulations to produce eruptions ([Linker et al.l 120031 : lAmari et al.ll2003l ). 



As it is observed on the Sun, flux cancellation is very ambiguous as it may be caused by sub- 



- 6 - 



mergence of fi-loops, emergence of U-loops, or magnetic cancellation by some dissipative process. 
Our simulation provides a unique opportunity to fully understand an example of flux cancellation 
in a way that is not possible through observations alone. We find flux cancellation occurs spon- 
taneously between opposite polarities that are driven together by the convective flows near the 
photosphere. The temporal evolution of cancellation is shown in Figure fTalfTc] which show a zoom- 
in view of the relevant area, with background color showing photospheric B z field. The horizontal 
motions in positive and negative polarities are shown by blue and red arrows, respectively. The flux 
cancellation event starts at time t = 05:00:00, lasts for 0.5 hr, and the total amount of cancelled 
flux approaches 1 x 10 20 Mx, which is 10% of the total unsigned flux on the photosphere. It is 
remarkable that during the process of flux cancellation, the coronal free energy (shown in the left 
panel of Figure [2]) is still increasing even though the photospheric magnetic flux is decreasing. 




Fig. 3. — : The structure of the magnetic fields in the area of flux cancellation at time t = 5:05:00 
from top (a and c) and side (b) view. The background planes show the photospheric B z field. Color 
on the rods indicates the vertical velocity at which the field lines are moving, with blue showing 
rising lines and red submerging. The red arrows in Panel (c) represent the horizontal magnetic 
field, with yellow line showing PIL. 



To understand this process, we take a close examination of the velocity fields in the area of 
flux cancellation prior to, during and in the later phase of the event, as shown by Figure fTa ] llbl and 
[Tel respectively. In Figure [Ta] the two polarities are pushed together by the horizontal converging 
motion at the two sides of the PIL. The flow pattern at the PIL is mostly converging and forms a 
very narrow PIL at a later time illustrated in Figure llbl Figure llbl shows a clear shearing motion 
at the PIL, with two polarities running in opposite directions. The strong shear flow at the PIL 
is driven by the Lorentz force, and takes place during the flux cancellation as the magnetic field 
expands in the corona as shown by the blue colored field lines in Figure [3a] and I3bl Figure da] shows 
the horizontal Lorentz force at time t = 05:01:00 with arrows, from which we find a clear pattern of 
the force running in the opposite directions across the PIL in the area of flux cancellation. At time 
t = 5:22:00, toward the end of the flux cancellation, the shearing motion is weakened and the flow 
pattern becomes dominated by the converging motion again, as shown in Figure [13 T he weak 



strong- weak shearing motion in the flux cancellation event is also shown by observations (ISu et al 
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Fig. 4. — : The structure of horizontal Lorentz force on Z = (a) and Z = 3 Mm (b) planes at times 
t = 05:01:00 (a) and 05:05:00 (b), respectively. The blue and red arrows represent Lorentz forces 
in positive and negative polarities, respectively. The magnitude of B z field component is shown in 
gray scale. Note that the Lorentz force reverses direction across the polarity inversion line driving 
the shear flow, and forms a spiral pattern in the negative polarity driving a rotating flow. 



20071 ) during flares where footprints of the magnetic fields present the same pattern of motion. 



Figure [3a] - [3c] show the configuration of the magnetic fields in the area of flux cancellation at 
time t = 05:05:00. Figure l3al provides the top view of the fields, and we find that instead of short 
loops perpendicular to the PIL, as is the case for potential fields, the field lines are instead elongated 
parallel to the PIL. Furthermore, the fields are more parallel to the PIL as they become closer to it, 
as shown by the red arrows in Figure [3c] At the photosphere, the magnetic fields are submerging to 
the convection zone at the PIL, as shown by the red rods in Figure [3b] The concentration of both 
the magnetic and velocity shear along the PIL demonstrates that the elongation of the field lines 
i s the result of shearing motion, which is further enhanced by the converging motion at the PIL 
dMartens fc Zwaan|[200ll ). These converging motions also gives rise to occurrence of t ether-cutting 



reconn ection close to the photosphere in the highly sheared fiel ds as proposed by 



Moore et al 



(|200ll ) as an explanation for solar eruptive events. Figure 1 of iMoore et al.l (120011 ) envisioned 
reconnection inside a highly sheared core field, which would produce longer, nearly horizontal field 
lines higher in the corona, while at the base of the arcade, a system of short unsheared loops 
would form. This flux transfer process accumulates the free energy of the sheared field higher in 
the corona while at the same time reducing the magnetic tension by detaching the system from 
the photosphere. The rise of the arcade naturally leads to necking off of the fields, which in turn, 
promotes more reconnection. Unabated, the runaway reconnection provides a mechanism for the 
initialization and growth of explosive events. With our simulation, we can fully examine the way 
in which submergence and reconnection alter the field line geometry in the tether-cutting process. 
A three-dimensional (3-D) view of the magnetic fields combined with the flow pattern provides a 
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complete picture of the effects of the internal reconnection. Figure I3bl illustrates the structure of 
the magnetic field lines colored by the vertical component of the velocity orthogonal to the field 
lines. Blue indicates upflows and red downflows. We find two groups of magnetic field lines formed 
during tether-cutting reconnection, with one rising up into the corona and the other, shorter loops, 
submerging into the convection zone. At the photosphere, the submergence of the shorter loops, 
as shown by the red rods right above the PIL, results in the decreased unsigned flux observed at 
the surface. Along the long, rising loops, magnetic shear accumulates, forming a highly sheared 
arcade. 
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Fig. 5. — : The structure of B 2 (a and d), F vert i C ai (b and e) and F emerg e (c and f) fields at the 
photosphere (a, b and c) and at Z = 3 Mm in the corona (d, e, and f) at time t = 05:05:00. PIL is 
shown by the black line. Red arrows in Panel (a) and (d) represent the horizontal magnetic fields, 
and white arrows in (b, c, e and f) show the horizontal velocity fields. 

To study the energy transport into the corona during flux cancellation, we examine the Poynt- 
ing flux at Z = 3 Mm associated with the vertical motions. At the surface, the vertical motion 
transports the energy back into the convection zone due to the concentration of the magnetic flux 
in downflows. In the corona, the acceleration of magnetic fields after reconnection plays a very 
important role in the energy flux, shown by Figure [5el At the PIL in the corona, we observe a 
strong energy input associated with the rising motion of the fields. The rising motion here in the 
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corona, however, is different from the emerging motion at the photosphere. At the photosphere, the 
emergence is caused by the expansion and buoyancy of the flux rope as well as convective upflows. 
However, in the corona, the tether-cutting reconnection forms longer loops possessing less magnetic 
tension, resulting in a magnetic pressure gradient force accelerating the magnetic fields upward. 




Fig. 6. — : 3-D structure of the magnetic fields in corona at time t = 5:05:00. The background 
planes show the B 2 field at Z = 3 Mm. Green rods represent low-lying field lines below 5 Mm with 
purple showing high field lines. The blue and red arrows show the horizontal velocity in positive and 
negative magnetic polarities, respectively. There is a clear shear flow along the polarity inversion 
line. The negative polarity also shows some signs of rotation, while the positive polarity shows only 
a coherent shear flow. The magnetic field lines clearly reflect the flow pattern. 

Besides reconnection and shear flows, rotation of the magnetic polarities also contributes to 
the accumulation of magnetic shear by moving the footpoints in opposite directions at the PIL. 
Figure [6] shows the geometry of two groups of coronal magnetic field lines: purple overlying lines 
with apex higher than 5 Mm, and green lines lower than 5 Mm. Horizontal velocity fields at Z = 3 
Mm, represented by the arrows, show a rotation pattern on the negative polarity, which is also seen 
in Figure O The rotation complements the shearing motion at the PIL, in the sense that both 
shear and rotation contribute to plasma motions parallel to the PIL. However the flow patterns 
are also quite distinct occurring at different scales and at different times. The shear flows are 
localized to the PIL and persist even when there is not a coherent rotation pattern in the magnetic 
polarities. When rotation is present, the shear flow tends to be greater in magnitude and will be 
shown graphically to dominate the energy transport into the corona. 
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The combination of the two motions creates a highly sheared arcade at the PIL, shown by 
the green rods in Figure [6j and rotation produces the twisted field structure at the far ends of the 
emerged fields, shown by purple rods. Both motions are driven by Lorentz force, as shown by the 
arrows in Figure llbl At the PIL, the horizontal Lorentz force runs in the opposite direction, consis- 
tent with the directions of the shearing motion. In the negative polarity, the Lorentz force rotates 
in the same way of the rotation of the polarity. The sigmoidal structure at the PIL is maintained by 
the consistent shearing and rotating motions in our simulation. The relative contribution of shear 
flows compared to rotation can be seen in the geometry of the coronal fields at the PIL. Here, the 
field is dominated by an arcade structure instead of a twisted flux rope. Most of the twist is located 
at the outer periphery of the two polarities. Even in the most twisted field lines, there is less than 
a full turn of the field around the axis of either polarity. The absence of the twisted flux rope may 
be explained by two facts: the axis of the initial magnetic flux rope remains in the convection zone; 
and there are not enough twisting motions on the surface to reproduce such a twisted structure 
after emergence. The emerged fields are sheared and elongated, forming the arcade structure over 
the PIL. 



4. Non-potentiality of the Magnetic Field 

During our simulation, a significant amount of energy, up to 2xl0 31 ergs, is transferred into 
the corona, over 5 hours of evolution. The question is then, how much of the magnetic energy is in a 
free form, which is necessary for solar eruptive events. To investigate the evolution of magnetic free 
energy, we compare the coronal magnetic fields in our simulation with potential fields extrapolated 
from the photospheric boundary conditions. Figure [7J shows the temporal evolution of the model 
fields (blue) and the extrapolated potential fields (red). Comparison between the two fields clearly 
shows the buildup of magnetic shear along the field lines during the process of flux emergence. At 
time t = 04:05:00, shown by the upper-left panel in Figure the magnetic flux concentrates in two 
polarities, with most of the field lines running perpendicular to the PIL. The energy transfer into 
the corona occurs as the magnetic field lines become elongated and sheared along the PIL. The 
lower-right panel in Figure [7] shows the field structure at time t =05:05:00, in which the field lines are 
sheared along the PIL and compressed into the lower atmosphere while the overlying fields remain 
almost potential. The simulated configuration of the magnetic fields in our doma in is consistent 



with observations and simulations of the magnetic fields before solar eruptions (jSchrijver et al 



2008 



Leake et al.ll2010i ). The overlying fields confine and compress the sheared core in the lower 



atmosphere, which may be destabilized and give rise to sudden release of magnetic free energy. 

Here, we calculate the free energy for the simulation by integrating the magnetic energy for 
the model and subtracting the energy of the potential field. The results are shown in the left panel 
of Figure [2] where the temporal evolution of model, potential and free energy is plotted with blue, 
red and black lines, respectively. Beginning at time t = 04:05:00, the amount of free energy doubles 
within one hour (the period shown by Figure [7j), such that the coronal free energy approaches 



Fig. 7. — : Comparison of the model field (blue) with the extrapolated potential field (red) at time 
t = 04:05:00, 04:25:00, 04:45:00 and 05:05:00. The color on the plane shows the B 2 field at the 
photosphere. The clearest departure in field line direction occurs at the sheared field lines above 
the PIL. 

7 x 10 30 ergs while the total energy is 1.7 x 10 31 ergs. At this time, the free energy contributes to 
40% of the total energy. Meanwhile, the vertical stratification of the free magnetic energy, is shown 
in the right panel of Figure [2j The line plot shows a strong tendency of concentrating free energy 
in the lower atmosphere, with 50% located below Z = 4 Mm, shown by the dashed line. Figure [Hal 
and I8bl shows the angle between the current and the magnetic field and the structure of a (the ratio 
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Fig. 8. — : (a) the angle between the current and magnetic fields (b) the a structure on the Y = 
plane. Black lines are the magnetic field lines and white regions in Panel (b) are areas where the 
angle between the current and magnetic fields are greater than 10°. 



of 1 4-7rj | and |B| when the angle between them is less than 10°) in the Y = plane, respectively. 
The non-force- free magnetic fields, as shown in Figure [5a| mainly reside in the convection zone, the 
core sheared field close to the surface and the negative polarity where the rotation forms a twisted 
magnetic structure (see Figure [7j). Note that a tends to be maximized low down the center of the 
arcade structure and then falls off with distance from the PIL. 



5. Poynting Fluxes 

To separate the energy transfer associated with vertical and horizontal flows, we decompose 
the Poynting flux into two components: 

-^horizontal = ~ "p - (B X U X + B y U y ) B z , (1) 
47T 

^vertical = ^{B 2 x +Bl) Uz . (2) 

Fvertical indicates the energy transported by the vertical motion, while Fhorizontai represents the en- 
ergy flux due to the horizontal motions, which include, at the photosphere, the separation of small 
dipoles, rotation, and the shearing flow along the PIL. Figure [5b] and [5c] show the structures of 
Fvertical and Fhorizontab respectively, with the arrows representing the horizontal velocity fields. It is 
clear that in regions with rotation or shearing there is a strong energy input into the corona. Com- 
parison of the structures of F ver tj ca j and Fhorizontai suggests that horizontal motion dominates in the 
energy transfer during the current phase of emergence, while vertical motion transfers energy back 



into the convection zone, due to the concentration of magnetic flux in downflow regions ([Fang et al 
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201dl . 12012 ). Vertical emerging motion dominates at the very begi nning of the flux eme r gence , and 
later gives way to horiz ontal motions, a pattern found earlier by iMagara &: Longcopd (|2003l ) and 
Manchester et al.1 (|200J). 



At the photosphere, the magnetic field is very highly structured by convection. In particular, 
the flow patterns and Poynting fluxes can be difficult to interpret with a large component of 
their structure being caused by random convective motions. To get a better understanding of the 
dynamics, we plot the vertical field strength, flow field and Poynting fluxes in the low corona at 
a height of Z = 3 Mm in Figure [5d] - [5fl Here, in the corona, the magnetic field expands and 
smooths out as the plasma beta drops. This transition allows us to make a much clearer picture 
of the relevant structures. The flux expands to form two nearby polarities and an elongated PIL 
between them, shown by Figure I5dl Here, the presence of a highly sheared field is obvious in the 
red arrows showing the horizontal field direction. The positive polarity shows little discernible twist 
while the negative polarity is clearly twisted. Flows are completely dominated by shearing motion 
at the PIL and rotation in the polarities. Both flows make a significant contribution to the energy 
transfer, but the energy transported by shear flow tends to be more consistent in time than that 
of rotational flows. 



6. Summary and Conclusions 

Our simulation provides us with a unique opportunity to study an example of magnetic flux 
emergence in a realistic simulation of the convection zone, coronal system. The system shows a 
wealth of complexity and interaction of mutiple physical processes, which conspire to transport 
magnetic flux and free energy from the convection zone into the corona. The emerging magnetic 
flux interacts with convective cells of varying scales as it approaches the photosphere, and then 
accumulates in convective downdrafts to form a bipolar magnetic structure resembling a solar active 
region. The magnetic field continually expands into the upper atmosphere, which results in shear 
flows and rotating motions that combine to draw the field nearly parallel to the PIL, and transport 
energy into the corona. At a later time, converging motions at the PIL cause flux cancellation at 
the photosphere along with tether-cutting reconnection, which produces highly sheared, sigmoidal 
shaped field lines high in the corona. 

In the process of building up the magnetic free energy, shearing flows play a significant role at 
the photosphere and in the corona. At the near surface layers, the angle between the current and 
magnetic fields ranges from 35° to 90° . The geometry of current and magnetic fields then produces 
the Lorentz force at the surface, which drives the rotating and shearing motions and increases the 
non-potentiality of the magnetic fields over a period of hours after the emergence. In particular, the 
shearing flows align the magnetic fields nearly parallel to the PIL and transport significant amount 
of magnetic ene rgy from convec t ion zo ne into the corona during the flux emergence, consistent with 



observations by lSchrijver et al.l (|2005l ). The shearing flows at the PIL are accompanied by periods 



of rotating motion of the magnetic polarities. Both of the two motions contribute to the build 
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up of magnetic shear and free energy in the corona. However, the shearing motions distinguish 
themselves from the rotation by their persistent presence and concentration at the PIL where the 
field has an arcade structure, while rotation on the other hand is only present in one of the magnetic 
polarities and lasts only for 1 hour. It is the long-lasting, Lorentz-force driven shearing motion that 
dominates the energy transfer. 

Convection-driven converging motion drives flux bundles of opposite polarities toge ther, pro- 
ducing strong magnetic g radient and the pronounced PIL, preferential for large flares (jSchriiver 



2007 



Falconer et al.ll2008l ) . The horizontal converging motion at the PIL also co ntributes to the de 



yelopment of the magnetic shear of the field lines by increasing the shear angle ([Martens k Zwaan 
200ll ). The combination of converging and shearing flows at the PIL forms a group of highly 
sheared arcades overlying a sharp magnetic PIL. Within this arcade, converging motions lead to 
the occurrence of tether-cutting reconnection, producing two types of field lines: one long sheared 
expanding loops and the other unsheared submerging loops. The short loo ps sink into the con - 
vection zone at the speed of convective flows, up to 2 km/s, consistent with lHarvey et al.1 ([1999), 
which continuously reduces the photospheric unsigned flux. The longer loops rise into the corona 
with more magnetic shear. The shearing motion along the PIL plays a very important role in 
that it both produces a sheared arcade structure ready for reconnection and accumulates magnetic 
shear in the field lines formed after the reconnection. The magnetic configuration of this area thus 
yields a high free ener gy up to 40% of th e total, which is c ompa rable with the threshold value for 
eruptions reported bv lAmari et al.l (|2003l ). lAulanier et al.l (|2010i ) and lMoore et al.1 (J2012J). 



Shear flows, converging motions, and tether-cutting reconnection combine to continuously build 
up the magnetic shear and free energy in the corona necessary for eruptive and explosive events. 
The magnetic reconnection and shearing motion at the photosphere produce the sigmoid-shaped 
field geometry, which is preferential for CMEs ( Canfield et al. Il999 ), and t he persistent sheared 



arcade structure is consistent with many CME models (e. g. 



Steinolfson 



1991 



Mikic k Linker 1994 



Amari et al. 



1996 



Guo k Wu 



1998 



Antiochos et al.lll999l ). Futhermore, we find coronal free energy 
grows at a rate of 3 x 10 30 ergs/hr, which b uilds up 10 31 ergs free energy over 5 hours. Such a 



build up rate is observed by ISchrijverl (|2007l ) prior to large flares. Moreover, the majority of the 
free energy resides low in the corona where it can be confined and later released by reconnection. 
Finally, we note that the consistent shearing motion and reconnection at the flux cancellation site 
keep reforming the sigmoid structure, which is essential for homologous eruptive events. 

In our simulation here, tether-cutting reconnection and flux cancellation take place in the 
highly sheared magnetic fields along the PIL, in the presence of converging and shearing flows. In 
the coupled system of the convection zone, the photosphere and the corona, all these mechanisms 
combine to work simultaneously and build up the free magnetic energy in the coronal fields. More 
importantly, the resulting geometry of the magnetic fields in our domain consists of compact, highly 
sheared core fields confined by more relaxed, overlying fields. 50% of the corona free magnetic 
energy is stored within the compressed, sheared core, while the overlying fields relax very quickly 
with altitude. The compact core, possessing free magnetic energy up to 10 30 ergs, provides the 
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preferential magnetic structure as well as free energy for solar eruptions. We expect that in future 
simulations of larger scales, shear flows, converging motion and reconnection will produce enough 
free magnetic energy and give rise to eruptions of the magnetic system. 
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